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Relationships serving as the basis for the calculation of optimum glass-annealing regimes are analyzed. The 
method for the calculation of optimum cooling regimes is extended to the general case of heat treatment of 
glass, including the heating, exposure, and cooling stages. 


The problems of identifying an optimum regime for glass 
annealing can be set in two ways. The first consists in mini¬ 
mizing the time the glass article stays in the annealing fur¬ 
nace under limited residual stresses. Such approach is justified 
in designing new furnaces, which may have a smaller length. 

Another setting of the optimum control problem aims for 
the minimum possible residual stresses in the glass articles 
under a preset duration of annealing in existing furnaces. 
However, the important purpose in technology is not to en¬ 
sure the minimum possible stresses but rather to keep stresses 
below the admissible values prescribed by regulations 
(50 - 100 nm/cm for rolled sheet glass and 50 - 200 nm/cm 
for glass containers). In the general case the annealing pro¬ 
cess, besides cooling, also includes heating and exposure of 
glass articles at a preset temperature. 

It is important to identify an annealing temperature 
schedule involving a minimum consumption of energy. This 
problem can be set as follows: for a preset duration of an¬ 
nealing At it is necessary to identify the dependence of the 
cooling rate on the temperature of the glass article R ( T ), un¬ 
der which a minimum temperature in achieved in the glass at 
the end of heating T w , provided that residuals stresses in the 
glass article do not exceed the permissible level. 

The author in [1] obtained a dependence of the cooling 
rate on the temperature of the glass article that ensures the 
minimum residual stresses in the article: 
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where A = J -yJW(T) dT; 7j and T h are temperatures signifi- 
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cantly lower and significantly higher, respectively, than the 
annealing interval. 

The temperature of a glass article can be taken as the 
temperature of a certain surface point (a glass band surface 
or an angle point on the bottom surface of a cylindrical arti¬ 
cle). The function W{T) characterizes relaxation processes 
in the glass and can be taken as the second derivative of the 
structural temperature. 

The structural temperature of a glass-forming material in 
an equilibrium state (T> 7j) is equal to the actual tempera¬ 
ture and in an isostructural state (T<7j) is constant [2]; 
therefore, the function W(T) outside the annealing interval is 
equal to zero, which according to expression (1) means infi¬ 
nitely high cooling rates. Consequently, an ideal annealing 
regime is impossible in practice but can be approached by 
cooling a glass article outside the annealing interval at the 
maximum rates possible for the particular furnace. 

All the above formulas and expressions are true for what 
is known as annealing “from above,” when an article enters 
the furnace with a sufficiently high temperature, under which 
all emerging stresses relax. In a real process glass often en¬ 
ters the furnace at an insufficiently high temperature corre¬ 
sponding to a nonequilibrium state of the glass-forming ma¬ 
terial, which prevents it from immediate cooling. Therefore, 
one should consider the general case of identifying an opti¬ 
mum annealing schedule, which apart from cooling includes 
stages of heating and exposure of glass articles. One should 
take into account the restrictions imposed on maximum heat¬ 
ing and cooling rates by the properties of the heaters and the 
lining of the annealing furnace. 

For a general case of an optimum annealing schedule the 
temperature history of a glass article in an annealing furnace 
should obey the following regularities. First, the glass has to 
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Fig. 1. Optimum annealing regime for rolled sheet glass: S. and S t 
are the residual stresses on the bottom and top surfaces of the glass 
band, respectively; S is the residual stress in the middles of the 
glass band; S is the maximum temporary stress; P h is the equiva¬ 
lent power of annealing furnace heaters. 


Otherwise it is impossible to implement the optimum an¬ 
nealing schedule. 

It is possible to use the “golden section” method in 
searching for the minimum temperature to end the heating of 
the glass article. The left boundary for this search is T- and 
the right boundary is a temperature obviously exceeding the 
vitrification temperature T h . This method converges to the 
sought-for r an in 7 to 10 steps. Simulation results revealed 
that the exposure stage is not necessary, i.e., after the article 
is heated to T an , it can be immediately cooled according to 
the optimum curve R(T). 

In calculating the optimum cooling curve, it is conve¬ 
nient to make all values discrete, i.e., to determine the values 
W( Tt) and7?(r ( .) for; = 1, m; (R(T t ) is the glass cooling rate 
in the temperature interval from T i to T i+] ; m = T an - T { ; 
T \ = r an > T i = r an “ l ’ •••> T m = T f + !)• The integral A in for¬ 
mula (1) in this case can be calculated using the trapezoidal 
rule: 

a) 
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be heated from its temperature at the furnace entrance T- up 
to the annealing temperature T with the maximum possible 
rate for this furnace A’. . The heating stage duration can 

“ max 

be found from the formula 

T - T 
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Then the article is exposed at the temperature T for the 
time ? exp required for the relaxation of stresses. The exposure 
stage ends when temporary stresses stop changing. After that 
the glass has to be cooled from its annealing temperature to 
the temperature of exit from the furnace T f based on a curve 
determined by the optimum cooling rates R(T) obtained 
from expression (1) taking into account the restrictions on 
the maximum cooling rate: 

R(T)<R r . 

max 


In calculating the optimum cooling rates R (T j ) based on 
expression (1), their values may exceed the maximum rate 
R that can be implemented in a particular annealing fur- 

nace. In this case one should take R (T,) = R r and update 
the values of the weight function W(T j ) using the formula 
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The integral of the root of the weight function A is calcu¬ 
lated from Eq. (2). Variation of the cooling rates will make 
the cooling time calculated for discrete velocities 
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For a furnace length / f and velocity of the article inside 
this furnace v art the total duration of annealing will be 

t - lf 
‘tot > 

'‘art 

and the cooling duration will be equal to 


For the average cooling rate the following relationship 
has to be satisfied: 

R m = Tm ~ T{ ' < R c . 

rn ^ ^ max 


different from the required At. Accordingly, the values have 
to be corrected until the difference between the absolute va¬ 
lues of f cool and A? becomes less than a certain value of accu¬ 
racy. 

As an example, let us calculate the optimum annealing 
schedule for rolled sheet glass at the Krasnyi Mai Glass 
Works (Vyshnii Volochok). The regime that is optimum for 
energy saving is shown in Fig. 1. 

The optimum regime was calculated using the automated 
system for calculation of glass-annealing schedules [3] for 
the following initial conditions: 7’ i = 450°C, T f = 250°C, 
l r = 61.8 m, v. = 1.2m/min; R„ = 15 K/min; 

R r =10 K/min. 

max 
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Using the optimum schedule allows for 8 - 10% saving 
in energy consumed in glass band annealing while maintain¬ 
ing the required quality level. The equivalent power of the 
heaters in this case is reduced due to a shortening of the glass 
exposure, whereas smoother cooling in the critical annealing 
zone facilitates a decrease in residual stresses. 
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